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Parity measurement is a key step in many entanglement generation and quantum error correction 
schemes. We propose a protocol for non-destructive parity measurement of two remote qubits, i.e., 
macroscopically separated qubits with no direct interaction. The qubits are instead dispersively 
coupled to separate resonators that radiate to shared photodetectors. The scheme is deterministic 
in the sense that there is no fundamental upper bound on the success probability. In contrast to 
previous proposals, our protocol addresses the scenario where number resolving photodetectors are 
available but the qubit-resonator coupling is time-independent and only dispersive. 
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I. INTRODUCTION 

One of the main challenges in quantum computing is 
that generally controlled qubit-qubit interactions should 
be strong while uncontrolled interactions should be neg¬ 
ligible. Gonceptually, the most straightforward approach 
to solving this issue is to physically remove uncontrolled 
degrees of freedom within a distance comparable to the 
qubit spacing, e.g., by trapping ions in ultra-high vac¬ 
uum |T]. An alternative approach is to entangle the 
qubits with photons that act as flying ancilla qubits. This 
approach allows placing the stationary qubits in remote 
locations because the photonic ancilla qubits generally in¬ 
teract weakly with the environment. In order to entangle 
the stationary qubits, the ancilla qubits only need to in¬ 
terfere optically with each other at detector inputs nig. 
Post-selection or local qubit operations conditioned on 
the detector outputs can then entangle the stationary 
qubits |2HT3]. Such entanglement generation is at the 
heart of quantum repeaters and cluster-state models of 
quantum computing mui]. 

More generally, the photon-mediated approach allows 
performing a remote non-destructive parity measurement 
(RNPM) [TT] described by the measurement operators 
ni = Igg) {gg\ + |ee) (ee| and n_i = \ge) {ge\ -b \eg) {eg\ 
for the even and odd parity outcomes, respectively. Here 
\g) and |e) are orthogonal single-qubit states that de¬ 
fine the computational basis. Parity measurement is 
more general than entanglement generation in the sense 
that RNPM can generate entanglement without ancilla 
qubits, but protocols for preparing entangled states do 
not necessarily allow measuring parity without the use 
of ancilla qubits. Specifically, Bell pairs are generated 
from an initial product state (|(jf) -f |e)) 0 (|g) -b |e)) /2 
by measuring its parity, which provides just the right 
amount of information to produce maximally entangled 
pairs. However, we envision that the most useful appli¬ 
cation for ancilla-free RNPM may be in quantum error 
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correction, where multi-qubit parity measurements are 
central to a wide variety of stabilizer codes |TS]. Re¬ 
mote parity measurement in particular may prove useful 
in circumventing limitations that are inherent to codes 
that only use geometrically local parity checks m- 

We note that the main challenge in non-destructive 
parity measurement is that the protocol must preserve 
coherence of arbitrary superpositions within the parity 
subspaces. In contrast, entanglement generation proto¬ 
cols drive the system into a known state and may assume 
a fixed initial state. In this sense RNPM is similarly chal¬ 
lenging as applying a GPHASE gate to remote qubits 
[7], although the two operations are not interchangeable 
without ancilla qubits. 

We propose a protocol for RNPM that can be im¬ 
plemented in circuit quantum electrodynamics (cQED) 
[nHzn using standard and minimalistic resources, with 
the exception of number resolving photodetectors with 
high temporal resolution. Such detectors for itinerant mi¬ 
crowave photons have not been realized to date but are 
under active theoretical and experimental development 
|22H31| . In addition to the detectors, the requirements 
for our protocol are the following: a beam splitter, two 
qubit-resonator systems with dispersive coupling, coher¬ 
ent drive pulses applied to the resonators, and single¬ 
qubit phase gates conditioned on the recorded arrival 
times of the photons. 

We emphasize that our protocol places few restrictions 
on the systems used as qubits. Specifically, we do not 
require A-type internal level structure to entangle the 
qubit state and a Fock state of the resonator, which 
is the typical starting point of most proposals for op¬ 
tical experiments [SHU El EH US] . Instead, our proposal 
works with a time-independent dispersive shift x ^-i^d res¬ 
onator decay rate k, both of which are standard features 
of cQED setups. The time-independent and dispersive 
qubit-resonator coupling gradually entangles the qubit 
state with coherent states of the resonator and contin¬ 
ues to play an important role throughout the protocol. 
This is in contrast to the generation of single photons on 
demand in Ref. [7], where the qubit state is entangled 
with Fock states effectively instantaneously compared to 
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other time scales. However, if on-off modulation of x is 
possible, we can use it to speed up the parity measure¬ 
ment by turning off the interaction at a specific time toff 
that maximizes the qubit-resonator entanglement. If the 
strong dispersive limit [35] y ^ /c is also reached, the in¬ 
teraction becomes effectively instantaneous compared to 
other time scales (tofi <C In this particular sce¬ 

nario, the protocol after toff reproduces a specific case of 
the protocol proposed by Azuma et al. in Ref. HU. 

Kerckhoff et al. proposed another parity measurement 
protocol that relies on homodyne detection and sequen¬ 
tial reflection of a probe signal from two resonators cou¬ 
pled to three-level atoms [33] . Roch et al. experimentally 
demonstrated the generation of odd-parity Bell states 
in cQED using a similar sequential setup, but with dis¬ 
persive qubit-resonator coupling and three distinct out¬ 
comes [34]. The disadvantage of the latter scheme is that 
it distinguishes the even parity states from each other and 
therefore reveals too much to function as RNPM. Our 
protocol on the other hand measures both parities non- 
destructively in a non-sequential setup. This is possible 
because the photodetectors effectively erase the phase in¬ 
formation that would allow distinguishing qubit states of 
the same parity. Instead, the photodetectors reveal the 
stochastic relative phase acquired by the states, which is 
the fundamental backaction of dispersive measurements 
[3SH37]. Conditioning phase gates on the photodetec¬ 
tor outputs therefore allows undoing the measurement- 
induced dephasing within the parity subspaces, much like 
in the extensively studied case of a joint measurement of 
two qubits in a single resonator [381130] . 

The remainder of this article is organized as follows. 
Section [T^ reviews measurement-induced dephasing and 
discusses its reversal in a single-qubit scenario. Sec¬ 
tion [Hljintroduces the RNPM protocol and demonstrates 
its validity under ideal conditions. It also suggests an 
alternative variation of the protocol for the strong dis¬ 
persive limit and briefly discusses some practical hurdles 
to implementing the protocol. Section [TV] concludes the 
article. 


II. REVERSING MEASUREMENT-INDUCED 
DEPHASING OF A SINGLE QUBIT 


Measurement-induced dephasing and the possibility of 
reversing it has been extensively discussed in the context 
of dispersive qubit measurement using quadrature detec¬ 
tors [351133] . The principle of reversing the dephasing 
using a photodetector has also been described by Frisk 
Kockum et al. in Ref. [35]. Nevertheless, we begin by 
reviewing these concepts in a single-qubit scenario be¬ 
cause it maps one-to-one to the even-parity subspace of 
the two-qubit scenario discussed in Sec. |III| 

We consider the qubit-resonator system illustrated in 
Fig. [3 We assume that the coupling is dispersive so that 
the closed-system Hamiltonian in the Schrodinger picture 



FIG. 1. (Color online) (a) Schematic diagram of a qubit dis- 
persively coupled to a resonator that radiates into a photode¬ 
tector. At t = 0, a short (Td ^ external pulse 

drives the resonator mode (a) from vacuum to a coherent 
state la). The resonator then evolves into a superposition of 
two coherent states at a rate x due to a dispersive interaction 
that entangles the qubit (a^) and the phase of the coherent 
state. At a rate k, the resonator decays through a single port 
to an output mode monitored by a photodetector that encodes 
the arrival time ti of each photon in its output record N (t). 
(b) The same setup depicted using cQED components. The 
optical cavity with asymmetric mirrors is replaced by a mi¬ 
crowave transmission line resonator with asymmetric coupling 
capacitors. The free-space links are replaced by transmission 
lines. The circulator ensures that the photodetector causes 
no direct back-action on the system, even if the detector is 
not ideal. 


is well approximated by 

Hq_r = (Wr -b X^z) o''O, + 

+ h[e {t) -b £* (f) e“'‘a] , 

as given in Ref. [^. Here, az = |e) (e| — \g) (gj, d is 
the annihilation operator for photons in the resonator, 
Wr/STT is the resonator frequency, (Iia/ 27 r is the Lamb- 
shifted qubit frequency, x/tt is the dispersive frequency 
shift of the qubit per photon, and e (t) describes the 
amplitude of a classical drive of the resonator through 
a weakly-coupled port. However, we work in an inter¬ 
action picture where we transform the basis states by 
exp [i (wrO^d -b uJadzI^j t] and correspondingly redefine 

iLq_r = hx^zd^d + h[e (t) -b s* {t) a] . (1) 

We consider a drive that displaces the resonator from vac¬ 
uum |0) into a coherent state |a) = D (a) |0) at t = 0 but 
is otherwise off. This is approximately the case for a short 
but strong Gaussian pulse e (t) = iae~^ j 
where k, x ^ ^ Wr. We note that the average pho¬ 

ton number |a| in the initial state is typically limited by 
assumptions underlying the dispersive coupling approxi¬ 
mation [36] . 

We assume that the dominant interaction between the 
qubit-resonator system and its environment is a weak 
linear coupling between the resonator and a transmis¬ 
sion line with a photodetector at the other end. We as¬ 
sume that the photodetector emits negligible noise power 
to its input and that it is impedance matched, i.e., the 
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detector is operated in the scattering mode [U]. Given 
the rotating-wave, secular, and Born-Markov approxima¬ 
tions, the corresponding stochastic master equation for 
the reduced density operator pi (t) of the qubit-resonator 
system conditioned on the photodetector output is 


dpi = dN (t) Q [y^d] Pi 
i 


— dtH 
dtT) 


■ 2 
Kd 


— a' a 


Pi 


PI. 


( 2 ) 


as given in Ref. m- Here, i] is the efficiency of 
the photodetector, dN (t) G {0,1} encodes the pho¬ 
ton arrival times T = {ti\dN (ti) = 1}, and {dN (t)) = 
dttr [ryKa^ap/] gives the detection probability. The su¬ 
peroperator Ti [c] p = cp + pc^ — tr \cp + pc^] p describes 
the continuous evolution between detection events, 
Q[c\p = cpc^/iv [cpc^] —p accounts for the discrete jumps 
at T, and I? [c] p = cpc'^ — | (c’^cp-f pc^c) includes the 
effects of unmonitored decay. The subscript I empha¬ 
sizes that the solutions, called trajectories, depend on 
the stochastic photodetector output. In general, Eq. ([^ 
should include terms corresponding to other imperfectly 
monitored decay channels, such as D p/ describ¬ 

ing spontaneous relaxation of the qubit, but we assume 
that K is sufficiently large to neglect them. 


A. Unmonitored system: measurement-induced 
dephasing 


In the limit p —>■ 0, Eq. ([^ reduces to the deterministic 
master equation 


dtp = 



+ V [v^a] p. 


( 3 ) 
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FIG. 2. (Color online) Example trajectory (solid lines) and 
average evolution (dashed lines) of the qubit-resonator sys¬ 
tem (see Fig. [T|) for y = k, a = 1, and the qubit initially in 
(Ip) -I- |e)) (a) Expectation values that are independent 

of the photon arrival times T = {ti}: square root of the pho¬ 
ton number d^d, the initially excited quadrature (a^ -|- a) /2, 
and product of the conjugate quadrature i (a^ — a) /2 and the 
qubit operator —dz. Individually, — a) = {dz) = 0 for all 
t (not shown), (b) Expectation values of d^ and dy which 
show decoherence of the qubit in the computational basis. 
The solid points a,t t — show the state after applying a 

phase gate R {(f>) that reverses the measurement-induced de¬ 
phasing (f> (T) [see Eq. (8 []. The curves describe numerical 
solutions to Eqs. and Q. 


Cegit) = aeg{t) / {agit)\aeit)) and Ceg (t) = c*g^{t), 
where {ag (t) [ue (t) ) = exp — |a|^ e~'^* (l — 


where we have dropped the subscript / since the evolu¬ 
tion is independent of the photodetector output. Fur¬ 
thermore, in the scattering mode the detector type is ir¬ 
relevant for determining p, and hence Eq. ^ is the same 
as for quadrature detection |361 [45] . 

Gambetta et al. solved Eq. ^ for a general drive 
e (t) [33]. For the specific parameters in our work, the 
solution can be expressed in closed form. For t > 0, it is a 
superposition of two coherent resonator states entangled 
with the qubit state: 

P(^)= H (t) I*) 01 ® (0) (ay (0I - (4) 

where ag = (cte = describes 

the exponentially decaying coherent resonator state given 
that the qubit is in \g) (|e)) [see Fig. I^a)]. The di¬ 
agonal elements Cgg and Cee remain at their initial val¬ 
ues since there are no terms in Eq. (§ that flip Ip) 
and |e). The off-diagonal elements decay according to 


aeg {t) = Ceg (0) exp 


= Ceg (0) exp 


-2ix J ae{t')a*g{t')dt' 

„ 1 _ p(2ix-K)t' 

-l“l T^WJ' 


( 5 ) 


and Ceg (0) is the initial value of the off-diagonal element. 

Here |ceg (t)| quantifies the phase coherence remaining 
in the qubit-resonator system as a whole, while 2aeg {t) = 
{dx {^v (0) describes the lateral components of the 

qubit Bloch vector after tracing out the resonator [see 
Fig.gb)[. For t ^ K the two are equivalent because 
the resonator state approaches vacuum exponentially ir¬ 
respective of the qubit state. Furthermore, in the long 


time limit Ceg (t) jceg (0) « exp 


l«lV(i 


*k/2x) 


so 


the net effect of the process on the qubit consists of a 
coherent rotation by an angle 


?f>uncond = - |q:|^ / (2x/«: + k/2x) (6) 


together with a reduction of coherence by a factor 
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of exp 


|a|V(l + «V4x^) 


The latter is called 


measurement-induced dephasing because it is directly re¬ 
lated to the ability to determine the qubit state by mon¬ 
itoring the radiation leaking out of the system [MUST]. 


B. Reversing dephasing using a photodetector 


In the opposite limit of perfect photodetection (ry 
Eq. Q reduces to 

dpi = dN (t) Q [v^a] pi 


— dtT-L 


* R -C 


d)d 


Pi, 


1), 


(7) 


which does not include any unmonitored decay channels 
and therefore does not change the purity of the initial 
state 133 . Since the purity does not change and there 
are no terms that flip \g) and |e), the initial qubit state 
must be restored in the long-time limit by a phase gate 

R{cli) = e-^^/^\g) (5|+e*^/2|e) (e|, 

where (j) is determined by solving Eq. Q for a given 
dN (t). Below, we show that this stochastic phase is in¬ 
dependent of the initial qubit state and can be written 
in closed form as 


^{T)=2xY.U. (8) 

ti^T 

Equation (|^ can be intuitively understood by noting 
that the photons are all injected into the resonator at 
t = 0, and hence the time each detected photon inter¬ 
acts with the qubit is equal to ti- The contributions add 
linearly so the total accumulated dispersive phase shift 
between |e) and |g) is as shown below for a pure 

initial state. Frisk Kockum et al. used an alternative 
method of solving the problem by applying a polaron 
transformation that allows writing a stochastic master 
equation for a two-level system only jl^. Specifically, 
choosing ry = 1, emtmeas = id, and tmeas < for 

the parameters defined in Ref. |35] corresponds to the 
scenario we consider. Note, however, that the analysis 
in Ref. |31] excludes the coherent evolution due to the 
ac Stark shift 2xRe (agCr*). Compared to our results, 
this leads to an additional phase factor that approaches 
exp (—i^uncond) for t 3> K~^ (see Eq. (26) in Ref. |32])- 

For a pure initial state, Eq. Q for the density 
operator pj = \'ij}i) ('0/| is equivalent to a stochas¬ 
tic Schrodinger equation for an unnormalized state 

Wi{t)) = US]- For t > 0, 

this means solving the Schrodinger equation for a non- 
Hermitian Hamiltonian /i — iK/2) affi between de¬ 
tection events, while a detected photon is taken into ac¬ 
count by applying the jump operator d and renormaliz¬ 
ing the resulting state. A photon is detected whenever 
(' 0 / (t) 107 (t)) reaches a random number drawn uni¬ 
formly and independently from [0,1] for each detection 
event i. 


For a normalized initial state {qg \g) qe\e))\o), the 

unnormalized state before the first detection event is 

10 ' ( 0 ) = 

-k qe|e)|aefo‘x-'"/2)‘ 

If Ti > , a photon is detected at ti = 


1 -I- |a| ^ln(ri) 


—K ^ In 
ter applying a becomes 


and the normalized state af- 


qge'‘^*^\g)\ae 


(iX-K'/2)t 


)q^e 




)• 


Similarly, each subsequent detection event results in fac¬ 
tors of leading to the normalized state 


10/ (f)) = % exp 




ti<t 


\9) 


qe exp 


-iX X! 


ti<t 


ae 


\e) 


(*X-K/2)t 


ae 


(-R-K/2)t 


where ti £ T and no more jumps occur when < 
exp |a|^ , with to = 0. See Fig. j^b) for an 

example with one detection event. 

Evidently applying R[(l){T\ti < t)] to |0/(0) undoes 
the relative phase between the \g) and |e) terms above. 
In the long-time limit t ^ «;“^log|Q!|, the overlap 


(0/(0) R(<()) 0/(0) 


furthermore approaches unity 


and additional detection events become exponentially un¬ 
likely as the resonator returns to |0). Therefore we con¬ 
clude that the combination of photodetection and R (0) 
effectively erases all the information that leaked out of 
the qubit during the process. In the parity measurement 
protocol, the same principle is used to erase only the 
part of the information that would allow distinguishing 
between \gg) and |ee) if a different scattering-mode de¬ 
tection scheme were used. 

Note that the resonator and the qubit periodically 
return to a product state, i.e., ag (ik) = Q^e {i^ = 
(—1)^ g-KR /2 g {kTrx~^\k G Z+} [see Fig.j^a)]. 

At these times the dynamics can be stopped by a second 
displacement D [—ag (^fc)] that brings the resonator to 
|0) deterministically, thereby decoupling the duration of 
the measurement from in the x > k regime. We take 
advantage of this possibility in the parity measurement 
protocol proposed in the next section. 


III. REMOTE PARITY MEASUREMENT 

Figure shows the setup we propose for measuring 
the parity of two remote qubits. It consists of 

two identical instances of the dispersively coupled qubit- 
resonator system described in the previous section. The 
two resonator modes oi and 02 are driven to an identical 
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coherent state |a) at t = 0 but the radiation leaking out 
of them is not measured individually. Rather, a 50:50 
beam splitter is arranged such that two identical pho¬ 
todetectors monitor the sum and difference modes with 
lowering operators c± = (hi ± 02 ) /\/2. The stochastic 
master equation for this system is 


dpi = dN+ (t) Q [y^c+] Pi 
+ dN_ (t) Q [y/rfRc^] pi 


— dtT-L 
+ dtV 
+ dtV 


+ fab. + f ala_ 

y/(l - r/) KC+ PI 
y/(l - ?7) KC- pi, 


Pi 


(9) 


where r] is the photodetection efficiency, N± (t) is the 
number of photons registered by the detector monitoring 
c±, and 


( 10 ) 

for t > 0. We denote coherent states of the sum and 
difference modes by |/3)j_. 

The parity measurement and the measurement- 
induced dephasing in this scenario are closely connected 
to the case of multiple qubits coupled to the same res¬ 
onator mode |38lH01 021 H5H5D] . Within the even-parity 
subspace, the difference mode in fact remains in vacuum 
so, mathematically, the description of the system maps 
one-to-one to the single resonator case. However, physi¬ 
cally the situation is distinct because the populated res¬ 
onator mode is the non-local sum mode. In the odd- 
parity subspace neither mode remains in vacuum. 


A. Protocol 

In this section, we describe the proposed RNPM pro¬ 
tocol and its effect on the qubits. Section [IlI B| proves the 
validity of these claims analytically for 77 —1. Addition¬ 
ally, Appendix ^presents some numerical trajectories in 
non-ideal cases. 

The steps of the protocol are the following: 



FIG. 3. (Color online) (a) Schematic setup for measuring the 
parity (Tz,\dz ,2 of two qubits located in identical resonators 
(see also Fig. [^. The resonators are driven into a coherent 
state I a) at t = 0 after which their phase becomes entangled 
with the local qubit at a rate X- At a rate k, the resonators 
radiate into a 50:50 beam splitter with photodetectors mon¬ 
itoring the output modes. The relative phases are chosen 
such that the detectors monitor c± = (hi ± 02 ) /\/2. At a 
final time ti (integer multiple of nx~^) the resonators are 
displaced back to vacuum. Recording more than zero pho¬ 
tons for c_ [N- {tf) > 0] indicates odd parity with certainty 
while N- (tf) = 0 partially projects the qubits to the even- 
parity subspace. The final N- (tf) is sufficient for reversing 
measurement-induced dephasing within the odd-parity sub¬ 
space but the full time-resolved A+ (t) is needed in the even- 
parity case, (b) The same setup depicted using cQED compo¬ 
nents. The beam splitter is replaced by a 180° hybrid coupler. 


5. Compute the measurement induced phases 


4'+ = ‘^x 

tieT+ 

and (/)- = ttN- (tf), 
and apply the local qubit feedback operations 


F = R 


2 2 


<R 


0+ <t>- 
2 2 


( 12 ) 

(13) 

(14) 


Assuming r/ 1, the state of the qubits after these steps 
is 


1. Start with the resonators in |0) and the qubits in 
an arbitrary pure state 


l^q(0)>= dijiij}- 

ije{e,g} 


( 11 ) 


2. Apply D (a) to both resonators at t = 0. 

3. Wait until tf = where k G TL'^ is arbitrary. 

4. Apply D — (—1)^ to both resonators. 


I V'q (tf)) = f y/i^H-l + IV'q (0)) 


(15) 


where Pi = 1 — P_i. The outcome of the parity mea¬ 
surement is indicated by P_i = ('(/'q(tf) 


n_i tjjq (tf) 

and is determined by N_ (tf) and P+= {t,|diV+(t,) = 1} 
according to 


P-i = 


' 1 if tV_ (tf) > 0 


\Qge\ + keg I 


cosyt, 

\tieT+ 


otherwise. 
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Evidently a single c_ detection event leads to com¬ 
plete parity projection to the odd-parity subspace 
= ~1)) while c_|_ detections in general lead 
to an exponential suppression of the odd-parity compo¬ 
nents, as shown in Fig. Therefore the outcome of the 
protocol is not strictly speaking binary, but rather con¬ 
tinuous in the interval P_i € [0,1]. However, in the 
limit of many photodetection events, the protocol is well 
approximated by a projective parity measurement with 
P_i S {0,1}, up to exponentially small corrections given 
by the exact result. The expected number of detected 
photons is jal^ (1 — and therefore tunable. Fur¬ 

thermore, if P-i is not sufficiently close to zero or one at 
ti, steps [^ through]^ of the protocol can be repeated arbi¬ 
trarily many times at the expense of increased operation 
time. In that case, a single feedback operation F applied 
after the last repetition should account for all detected 
photons. The displacements in steps and of adjacent 
repetitions can be combined into a single operation as 
well. 

Finally, we claim that the outcome probabilities 
are distributed according to the parity of the initial 
state. More specifically, the ensemble average of P_i 


is E[P_i] 


, (0) 


n_i 


V'q (0) , 


This implies that 


the protocol is indeed a parity measurement, rather than 
some other operation that leads to a well-defined parity. 


B. Full temporal evolution 


Let us first give a qualitative explanation for Eqs. (121- 


(151 by rewriting the Hamiltonian in Eq. (101 in the basis 
of the monitored operators as 


i/(2) = ^ 
q-i- 9 


(CTz.I + (T^,2)(c+C+ -I- cLc-) 

+ (0-z,l - CTz,2)(c+C_ -f C+C'_) 


Note that after the initial displacements the sum and 
difference modes start in \^/2a)^ and |0)_, respectively. 
Furthermore, there are no terms in Eq. ([^ that flip the 
qubits in the computational basis so parity is conserved 
if it is initially well defined. For an even-parity initial 
state, (Tzq — 172^2 yields zero so there are no terms in 
Eq. (j^ that excite the difference mode out of |0)_ and 
we can therefore trace it out without loss of information. 
The remaining terms are identical to the single qubit case 
[see Eq. 01 with the mapping \gg) —>■ j^), jee) —> |e), and 
c+ —>■ a. This explains why Eq. (12) matches Eq. ([^. On 
the other hand for an odd-parity initial state, (Tz,i + <5’z,2 

t ^ (2) 

yields zero so the dispersive clj_c± terms in iLq.f vanish. 
Instead, photons are exchanged between the two bosonic 
modes by the clj_Cqi terms with a phase flip between \ge) 
and \eg) associated with each exchange. Since the only 
other event that changes the photon number in c_ is a 
detection event dN_ = 1, the parity of N_ (tf) and the 
parity of the number of phase flips must match given that 



t (units of K t (units of k 


FIG. 4. (Color online) Example trajectories (solid lines) and 
average evolution (dashed lines) for the two-qubit setup shown 
in Fig.[^with y = k, a = 1, tf = and the qubits initially 

in (|g) + |e))®2 /2. The left column (a,c) shows a trajectory 
that leads to projection of the qubits into the odd-parity sub¬ 
space, while the right column (b,d) shows partial projection 
into the even-parity subspace. (a,b) The expected photode¬ 
tection rates in units of k for the sum and difference modes 
(c±), with indicating the first and only photodetection 
event in the selected trajectories. Up to t{, the local field 


quadratures ^at -|- aj^ /2 and /2 are identi¬ 

cal to those shown in Fig.[^a), regardless of detection events. 
(c,d) Expectation values of the Bell state stabilizers 
where i G {x,y,z}. A projective parity measurement of the 
specified initial state leads to (dz, 1 ( 72 , 2 ) = ±1 and conserves 
{dx,idx, 2 ) = 1. The solid points at t = show the state 
after applying a phase gate F [see Eq. (141] that reverses the 
measurement-induced dephasing, i.e., restores ((jj, 1 ( 72 ;, 2 ) = 1. 
The curves describe numerical solutions to Eqs. 0 and (161 
with g = 1. Overlapping points have been offset horizontally 
for clarity. 


c_ starts and ends in vacuum and 77 = 1. This explains 

Eq. m. 


Let us formally prove the validity of Eqs. (12|-(15l 


and the associated probabilities for a pure initial state 
by explicitly solving the stochastic Schrodinger equation 
corresponding to Eq. ([^ in the 77 —>■ 1 limit. Analogously 
to the single-qubit case, the jump operators are c± and 
the non-Hermitian Hamiltonian that determines the evo¬ 
lution of the unnormalized state (t)) is 

H' = — ih'^ (^Cq_c+ -I- c[lc_^ 

i=i 


for t > 0. A photon is detected whenever {'tpj (t) {ipj (t)) 































7 


reaches a random number drawn uniformly and inde¬ 
pendently from [0,1] for each event i. The detector that 
clicks is chosen at random with probabilities weighted by 


For the initial state given in Eq. (Ill, the unnormalized 
state IV'j (t)) before the first detection event is 


+ qee\ee)\Ae~"^*) _^_\0) _ 

+ qge\ge)\Acosxt) ^\iAsmxt) _ 

+ <leg\eg) \ A cos xt)^ \ -iAsinxt)_'j , 


where ti G T+. At t( = kTTx~^, the difference mode re¬ 
turns to vacuum and the sum mode is in | (—1) A (t) 
The sum mode is then driven back to vacuum by the sym¬ 
metric displacements in step|^of the protocol. Therefore 
the dynamics stop and the resonators can be traced out 
without loss of information. 

In the N- (tf) > 0 case, the contribution of (/>_ in F 
[see Eq. (14l] undoes the relative factor be¬ 

tween \ge) and \eg), while the contribution due to (/)+ does 
nothing. Therefore IV'q {U)) — F [(0+0_ IV'/ (tf))] indeed 
matches Eq. (15) with P_i = 1. In the N_ (tf) = 0 case, 
the contribution of (j)^ in F restores the initial relative 
phase in the even-parity subspace so that 


where A (t) = y/2ae is an exponentially decay¬ 
ing amplitude of the remaining radiation. If ri > 


exp 


- 2 |a|^ (1 - 


, a photon is detected at ti = 


In 


l-k|a| ^ln(ri)/2 


the event is proportional to 


and the state \ipi (t)) after 


iqge I ge) IA cos ) _^ I * A sin xt) _ 

+ {-'i) qeg\eg)\Acosxt)^ \ -zAsinyt) 


if dN- (tf) = 1, or 


l^q {ti)) OC qgg\gg) + ^ee|sc) 

+ cosx^i^ {gge\ge) +qeg\eg)) . 

Note that in the limit iV_|_ (tf) 3> 1, the product of the 
cosines in general vanishes exponentially. 

To show that E[P_i] = (V'q (0) |n_i|'0q(O)) we use 
the fact that the expectation value of ll_i for the uncon¬ 
ditioned system state p is invariant in time. The time 
invariance follows from the fact that the master equation 


e^^^-qgg\gg)\AF^*)^\Q)_ 

+e-*^‘^gee|ee)|Ae-*^‘)^| 0 )_ 

-f cos (xti)gge|ffe)|Acosxf)^|*Asinxt)_ 
+ cos (xt 1) geg I e^) IA cos Xt)q_ | - sin xt) 


dtp =-^ >/3 [Vkc±] p 


I 

h L 


± 

2 




[CkOj] p 
i=i 


(16) 


if dN+ (tf) = 1. See Fig.|^for examples. The calculation 
for each subsequent event at tf is identical given the new 
initial state at tf_i. Therefore, at an arbitrary time t < 
tf, the state is 

\il)i (t)) OC t^-(*^gge|5e)|Acosxt)^|iAsinxt)_ 

+ (-t)^“^‘^ geg|eg)|Acosxt)_^| - tAsinxt)_ 
if N- (t) > 0, and otherwise 


IV’7 (t)) OC exp qgg\gg)\AF^^)^\t))_ 

-kexp 9 ee|ee)|Ae"*^*)^| 0 )_ 

+ cosxt*^ 

X ^gse|5e)|Acosxt)_^|tAsinxt)_ 

+ g.g|e,)|Acosxt)J-asinxt). 


does not depend on the type of scattering-mode detec¬ 
tion and is solved by p = piq ® pig, where p^ (t) is 
the solution to Eq. ([^ given in Eqs. (|^ and ([^. As 
noted in Sec. II A| the diagonal elements of piq (t) do 
not change in the computational basis. Therefore the 
diagonal elements of the two-qubit p {t) also remain con¬ 
stant, and hence tr p {t) n_i is time invariant. Ap¬ 
plying this time invariance to p (0) = ['i/’q (0)) (V'q (0)| 
shows that (i/jq (0) |n_i IV'q (0)) =tr p(0)n_i is equal 
to tr /3(tf)n_i =E[P_i]. 


C. Tunable coupling 

Here we propose a variation of the protocol for a sce¬ 
nario where the dispersive coupling can be effectively 
switched on and off in time, either by changing x or by us¬ 
ing dynamical decoupling m- If this capability is used to 
turn off the unitary evolution {h)^} —t 0) at toff = tt/Sx, 
the results of the previous section are valid with the re¬ 
placement of xt by X min (t, toff)- Consequently, a single 
detection event after toff leads to complete parity projec¬ 
tion because odd (even) qubit parity is associated with 






















the sum (difference) mode being in vacuum. In this mod¬ 
ified protocol we also require tf ^ k~^ and skip stepof 
the protocol. 

This variation of the protocol is particnlarly beneficial 
in the strong dispersive limit |35] x ^ ^ where most pho¬ 
tons leak out of the resonators after toff- In addition to 
complete parity projection, the variation may be of prac¬ 
tical benefit because (j)+ [see Eq. @1 becomes indepen¬ 
dent of the arrival times of the detected photons if they 
can be assumed to all arrive after toff, i.e., if x ^ \ct\^ k. 
In the extreme limit k/x 0, the initial time interval np 
to toff can be approximated as an instantaneons entan¬ 
glement operation between each qnbit and its resonator. 
In this special case, our protocol after tos coincides with 
the protocol proposed by Aznma et al. for the initial 
state 9 = TT, where 6 is defined in Ref. mi. 

A detection event at t > toff also projects the resonator 
into a known coherent state, so in principle the require¬ 
ment tf ^ can be relaxed by driving the resonators 
into vacunm by displacements conditioned on N±. Alter¬ 
natively, X can be restored to its initial valne at some time 
t' so that the difference mode evolves back to vacnnm at 
t' + toff regardless of the qubit state. At t' -|- toff uncon¬ 
ditioned displacements on c_|_ can then stop the process 
analogously to step of the original protocol. 


D. Practical considerations 


The main experimental hurdle to implementing onr 
protocol in cQED is that it reqnires nearly ideal pho¬ 
todetectors. Specifically, in order for the measurement to 
be non-destructive the photodetectors need to have high 
quantnm efficiency, low dark-connt rate, photon number 
resolution and, for the even-parity outcome, high tem¬ 
poral resolution compared to x~^- Here we use number 
resolution to mean that the detector must not have a 
significant dead time after a detection event. By high 
qnantnm efficiency we mean that, in addition to high 
detector efficiency, photon losses in other parts of the 
setup must be negligible. Failing to satisfy any of these 
requirements leads to erroneous terms in Eq. (12 \ and 


therefore randomizes the relative phase within the par¬ 
ity snbspaces. Appendix |A] shows the effect of imperfect 
qnantnm efficiency on example trajectories. 

It is possible to relax some of the above requirements: 
Temporal resolntion is nnnecessary in the x ^ ^ case if 
nsing the dynamical decoupling discnssed in Sec. [Ill C| 
Nnmber resolution becomes unnecessary if |a|^ <C 1 and 
the protocol is instead iterated many times. However, 
the latter increases the duration of the protocol and will 
eventnally invalidate the assnmption that other relax¬ 
ation mechanisms are negligible. For reference, x/^t'" in 
cQED is often several megahertz so a single iteration may 
take t{ ~ 100 ns. This should be compared to qubit co¬ 
herence times that have recently reached ronghly 100 /is 
[5^1 - 15^ . Appendix [A| presents some example trajectories 
for non-negligible qubit relaxation. 


The assumption of identical qubit-resonator systems is 
another source of concern for practical implementations. 
In cQED, X is nsnally tnnable throngh the qnbit fre- 
qnency bnt, typically, the resonator freqnency Wr and the 
decay rate k are not tunable. Fortnnately, the parameters 
of typical cQED resonators are highly reproducible |55| . 
Furthermore, in-situ tuning of both and k is possi¬ 
ble at the cost of increased complexity Finally, 

choosing an asymmetric drive in step [^ of the protocol 
can compensate for different values of n. In general, this 
only adjnsts the average nnmber of emitted photons per 
resonator and not the time scale of their emission. How¬ 
ever, in the strong dispersive limit it is possible to choose 
tf ^ and approximate the photon emission rate as 
constant. 


IV. CONCLUSION 


We proposed a protocol for remote non-destructive 
parity measurement of two qubits. The protocol is de¬ 
terministic in the sense that it leads to complete par¬ 
ity projection with a probability that approaches unity 
in the ideal case. Furthermore, it conserves the relative 
phase within the parity snbspaces even when the parity 
projection is incomplete. Therefore the protocol is also 
of repeat-until-snccess type [7] in the sense that it can 
be repeated until the desired degree of parity projection 
is reached. We proved these claims analytically for the 
ideal case and investigated effects of some of the practical 
limitations nnmerically. 

Except for requiring high-qnality photodetectors, onr 
protocol is experimentally implementable in cQED with 
minimalistic resonrces. In particular, the protocol places 
few requirements on the qubits and their control lines as 
it requires only time-independent and dispersive qubit- 
resonator conpling. This is promising for scalability and 
calls for a fnture extension of the protocol to many- 
qnbit scenarios. Such an extension would further reduce 
the overhead of measuring non-local mnlti-qnbit parity 
checks for the purposes of quantum error correction. 
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FIG. 5. (Color online) Trajectories with imperfect detection 
{rj = 0.9 for both detectors), which leads to imperfect reversal 
of the measurement-induced dephasing. See Fig. for the 
values of other simulation parameters and an explanation of 
the symbols. 



FIG. 6. (Color online) Trajectories in the presence of qubit 
relaxation at a rate 0.1 x k/tv. See Fig.j^for the values of other 
simulation parameters and an explanation of the symbols. 


Appendix A: Detector inefficiency and qubit 
relaxation 


Here we briefly discuss the effect of imperfect photode¬ 
tection and of qubit relaxation. We point out some is¬ 
sues that inevitably arise in experimental realizations but 
do not attempt to thoroughly map the effects of non¬ 
idealities. We do this by presenting numerical solutions 
to Eq. 


Figure presents example trajectories for imperfect 
detection efficiency {rj = 0.9). The possibility of 
missing photons leads to a mixed state and prevents 
fully reversing the measurement-induced dephasing, i.e., 
< 1 at f = ff. However, the unconditioned 
master equation [Eq. (16)] is unchanged and the parity 
of the initial state is correctly measured, as long as suffi¬ 
ciently many photons are detected. 


Figure shows trajectories for a non-zero qubit relax¬ 
ation rate. Specifically, we add T> [^/7i7_ j] pi for each 
qubit j to the right-hand side of Eq. ([^. We choose 
7 = 0.1 X k/tt. As in the case of inefficient detec¬ 
tion, the qubit relaxation prevents perfect reversal of the 
measurement-induced dephasing. In addition, qubit de¬ 
cay leads to a mixed resonator state even at tf. This 
implies that the displacements in step [^ of the protocol 
cannot restore the resonators to vacuum perfectly, i.e., 
(c!|_c±) > 0 even after t = t(. Visually the most striking 
phenomenon in Fig. is the non-conservation of parity 
but it occurs even without performing the measurement, 
i.e., even if we were to choose a = 0. 
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